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The  main  objective  of  this  paper  is  to  continue  the  development  of  activities 
of  basic  and  applied  research  related  to  wind  energy  and  to  develop  methods 
of  optimal  control  to  improve  the  performance  and  production  of  electrical 
energy  from  wind.  A  new  control  technique  of  Double  fed  induction 
generator  for  wind  turbine  is  undertaken  through  a  robust  approach  tagged 
nonlinear  sliding  mode  control  (SMC)  with  exponential  reaching  law  control 
(ERL).  The  SMC  with  ERL  proves  to  be  capable  of  reducing  the  system 
chattering  phenomenon  as  well  as  accelerating  the  approaching  process.  A 
nonlinear  case  numerical  simulation  test  is  employed  to  verify  the  superior 
performance  of  the  ERL  method  over  traditional  power  rate  reaching 
strategy.  Results  obtained  in  Matlab/Simulink  environment  show  that  the 
SMC  with  ERL  is  more  robust,  prove  excellent  performance  for  the  control 
unit  by  improving  power  quality  and  stability  of  wind  turbine. 
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1.  INTRODUCTION 

Wind  energy  is  a  free,  renewable  resource,  so  no  matter  how  much  is  used  today,  there  will  still  be 
the  same  supply  in  the  future.  Wind  energy  is  also  a  source  of  clean,  non-polluting,  electricity.  Unlike 
conventional  power  plants,  wind  plants  emit  no  air  pollutants  or  greenhouse  gases  [1], 

The  latest  generation  wind  turbines  operate  at  variable  speed.  This  type  of  operation  makes  it 
possible  to  increase  the  energy  efficiency,  to  lower  the  mechanical  loads  and  to  improve  the  quality  of  the 
electrical  energy  produced. 

A  Doubly  Fed  Induction  Generator  (DFIG)  is  an  electrical  asynchronous  three-phase  machine  with 
open  rotor  windings  which  can  be  fed  by  external  voltages.  The  typical  connection  scheme  of  this  machine  is 
reported  in  Figure  1 .  The  stator  windings  are  directly  connected  to  the  line  grid,  while  the  rotor  windings  are 
controlled  by  means  of  an  inverter  [2],  This  solution  is  very  attractive  for  all  the  applications  where  limited 
speed  variations  around  the  synchronous  velocity  are  present,  since  the  power  handled  by  the  converter  at 
rotor  side  will  be  a  small  fraction  (depending  on  the  slip)  of  the  overall  system  power.  In  particular,  for 
electric  energy  generation  applications,  it  is  important  to  note  that  the  asynchronous  nature  of  the  DFIG 
allows  producing  constant  frequency  electric  power  with  a  variable  mechanical  speed,  in  addition  reduced 
copper  losses  and  wider  operational  range  are  obtained  with  respect  to  standard  squirrel-cage  induction 
machine  [3], 
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The  sliding  mode  control  (SMC)  approach  has  heretofore  been  used  for  many  dynamic  control 
systems  due  to  its  simplicity  of  implementation  and  their  robustness  against  both  plant  uncertainties  and 
external  disturbances,  suitable  transient  performance  and  quick  response. 

The  SMC  has  gain  considerable  concerns  in  the  application  of  DFIG  control.  In  [4],  a  dynamical 
sliding  mode  power  control  scheme  is  proposed  for  DFIG.  In  [5],  the  SMC  method  is  employed  to  control  the 
rotational  speed  in  the  rotor  side  as  well  as  the  dc-link  voltage  in  the  grid  side  of  a  wind  energy  conversion 
system;  a  model  predictive  control  based  SMC  is  proposed  to  control  a  three-phase  grid-connected  converter 
and  the  grid  current  total  harmonic  distortion  and  the  switching  losses  are  largely  mitigated  in  [6]. 

One  major  bug  of  using  SMC  in  practical  applications  is  the  chattering.  Chattering  causes  several 
damages  to  the  mechanical  components.  Different  methods  were  proposed  to  minimize  chattering  such  as  by 
replacing  the  sign  function  with  saturation  [7]  or  using  Quasi-SMC  and  fuzzy  mathematics  [8].  In  [9]  a 
optimized  sliding  mode  control  (SMC)  strategy  to  maximize  existence  region  for  single-phase  dynamic 
voltage  restorers.  It  is  shown  analytically  that  there  exists  an  optimum  sliding  coefficient  which  enlarges  the 
existence  region  of  the  sliding  mode  to  its  maximum.  Also,  it  is  pointed  out  that  the  optimum  sliding 
coefficient  improves  the  dynamic  response.  In  [10]  the  author  proposes  a  SMC  control  strategy  which 
combines  a  switched  policy  with  a  time-based  adaptation  of  the  control  gain,  thereby  allowing  to  effectively 
deal  with  a  very  conservative  prior  knowledge  of  the  upper  bounds  on  the  uncertainties,  that  usually  leads  to 
a  large  control  authority,  and  related  performance  degradation.  In  [11]  this  paper  proposes  two  model -free 
sliding  mode  control  system  structures.  The  new  structures  are  compared  with  a  model-free  intelligent 
proportional-integral  (PI)  control  system  structure.  The  design  and  implementation  of  adaptive  Second  Order 
Sliding  Mode  Control  for  a  four  wheels  Skid-Steered  Mobile  Robot  [12]. 

The  combination  of  a  sliding  mode  control  algorithm  and  fuzzy  logic  developing  for  robot 
manipulator’s  trajectory  [13].  This  paper  aims  to  propose  a  SMC  control  scheme  based  on  a  novel 
exponential  reaching  law  (ERF),  the  objective  of  this  paper  is  the  control  of  active  and  reactive  powers  in  the 
variable  speed  wind  turbine.  The  nonlinear  sliding  mode  control  with  an  exponential  reaching  law  is  used  to 
achieve  our  objective.  In  [4]  we  has  used  the  conventional  sliding-mode  control  to  control  the  power 
produced  by  a  wind  turbine,  but  in  this  paper,  we  will  use  the  same  method  but  by  employing  an  exponential 
reaching  law  (ERF)  proposed  by  Fallaha  et  al.  to  calculate  the  control  law  [14]. 

The  advantage  of  the  proposed  approach  is  that  allows  to  obtain  a  good  reference  tracking  and  the 
suppression  of  chatter  phenomenon.  This  paper  is  organized  as  follows.  Firstly,  In  Section  2,  the  SMC 
conventional  and  SMC  with  ERF  methods  were  be  introduced.  In  Section  3,  the  ERF  based  SMC  will  be 
utilized  in  the  control  of  DFIG  for  wind  turbine  generator.  In  Section  4,  simulation  will  be  conducted  and  the 
performance  of  the  proposed  method  will  be  tested  with  the  simulation  results.  In  Section  6,  the  conclusion  of 
the  entire  paper  will  be  provided. 


2.  THE  CONTROL  SCHEME 

In  this  section,  first  we  represent  the  theoretical  structure  of  the  simple  sliding  mode  and  also  the 
exponential  sliding  mode  control  for  the  general  dynamics  of  rigid  systems. 

2.1.  Sliding  Mode  Control  (SMC) 

The  SMC  was  widely  used  in  the  literature.  This  success  is  due  to  its  simplicity  of  implementation 
and  to  its  robustness  towards  the  parameter  variations  and  the  external  disturbances.  It  is  to  be  switched  with 
the  discontinuous  functions  of  the  dynamic  system  the  structure  so  that  the  state  vector  to  follow  a  trajectory 
S(X)  =0  in  the  state  space. 
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2.1.1.  Choice  of  switching  surface 

The  surface  of  sliding  is  a  scalar  function  such  as  the  error  on  the  variable  to  adjust  glide  to  this 
surface  and  aims  towards  the  origin  of  the  plan  of  phase.  So,  the  surface  represents  the  wished  dynamic 
behavior.  We  find  in  the  literature  of  various  forms  of  the  surface,  among  which  each  gives  better 
performances  for  certain  uses.  J.  Slotine  proposes  a  shape  of  general  equation  to  determine  the  surface  of 
sliding  which  insures  the  convergence  of  a  variable  towards  its  wished  value,  [5]: 

S(X)  =  (^+Axy-le(X)  (1) 

at 


where: 

r  is  the  relative  degree  of  the  sliding  surface. 
e(x)  is  the  erreur 

Xx  is  a  constant. Concerning  the  development  of  the  control  law. 

The  parameters  of  the  sliding  surface  will  be  chosen  so  that  the  error  trajectories  will  start  from  the 
initial  point  to  the  sliding  surface  and  continue  to  slide  near  the  sliding  surface.  The  first  stage  is  called  the 
reaching  stage  and  the  second  stage  is  called  the  sliding  stage.  In  the  sliding  stage,  the  error  trajectories  will 
follow  the  following  Lyapunov  condition, 

V{X)  =  ^S2(X)  (2) 

Its  derivative  given  by: 

y(x)  =  s(x)s(x)  (3) 

2.1.2.  Convergence  condition 

One  must  decrease  of  the  Lyapunov  function  to  zero.  For  this  purpose,  it  is  sufficient  to  assure  that 
is  derivative  negative.  So: 

S(X)S(X)  <  -n\S(X)\  for  VX  *  0  (4) 

where  r)  is  a  strictly  positive  constant.  Satisfying  above  condition  keeps  the  system  trajectories  remaining 
on  the  sliding  surface. 

One  need  to  note  that  the  S  involved  in  equation  (4)  is  different  from  the  S  stated  in  equation  f  1). 
The  S  in  equation  (4)  stands  for  the  reaching  law  while  the  S  in  equation  (1)  represents  the  sliding  surface. 
Apparently,  it  can  be  seen  from  equation  (4)  that  S  and  its  derivative  S  have  different  signs  with  each  other, 
and  that  is  the  mean  cause  of  the  chattering  phenomenon.  Suppose  S  is  positive  at  the  starting  point,  then  its 
derivative  is  negative,  which  makes  S  approaches  the  sliding  surface  where  the  value  of  S  is  equal  to  0; 
when  S  reaches  the  sliding  surface  for  the  first  time,  it  will  continue  to  decrease  until  its  value  becomes 
negative  and  its  derivative  becomes  positive,  hence,  the  trajectory  will  increase  subsequently  until  it  becomes 
positive  again  and  similar  process  repeats  again.  That  is  the  cause  of  the  chattering  phenomenon. 

Two  parts  have  to  be  distinguished  in  the  control  design  procedure.  The  first  one  concern  the 
attractivity  of  the  state  trajectory  to  the  sliding  surface  and  the  second  represents  the  dynamic  response  of  the 
representative  point  in  sliding  mode.  This  latter  is  very  important  in  terms  of  application  of  non-linear  control 
techniques,  because  it  eliminates  the  uncertain  effect  of  the  model  and  external  perturbation.  Among  the 
literature,  one  can  chose  for  the  controller  the  following  expression: 

ll  —Ueq+Un  (5) 

where:  uec/  the  control  function  is  defined  by  Utkin  and  noted  equivalent  control,  for  which  the  trajectory 
response  remains  on  the  sliding  surface  [15],  un  is  the  discrete  control  which  is  determined  to  verify  the 
convergence  condition. 
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un  =  K  sigt{s(x ))  (6) 

Two  major  problems  exist  within  the  study  of  SMC  theory:  the  chattering  phenomenon  and  the 
approaching  speed,  which  largely  impact  the  application  of  SMC.  The  appearance  of  the  chattering 
phenomenon  near  the  sliding  manifold  is  inevitable  as  long  as  the  high  speed  switching  law  is  utilized.  Thus, 
the  appearance  of  chattering  phenomenon  restricts  its  domain  of  application  due  to  the  fact  it  may  excite 
some  undesirable  high  frequency  dynamics  which  further  causes  severe  damages  to  the  system.  Literature 
[16]  proposes  a  kind  of  ERL  based  SMC  technique  which  can  reduce  the  chattering  phenomenon  to  a  certain 
extent.  In  order  to  further  accelerate  the  reaching  speed  from  the  initial  point  to  the  sliding  manifold,  a  novel 
SMC  method  based  on  ERL  will  be  proposed  in  this  section. 

2.2.  Sliding  Mode  Control  with  Exponential  Reaching  Law 

With  the  aid  of  a  new  approach  of  choosing  K  we  can  control  the  tracking  error  as  well  as 
preventing  the  chattering  occurrence.  We  can  select  K  by  an  exponential  procedure  [17]. 

S=^—sig?7{S)  (7) 

N(S)  5 

where 

N(S)  =  j3+a-j3)e~r^  (8) 

a,  and  ju  are  all  positive  numbers,  |S|  denotes  the  absolute  value  of  5  ,  and  0  <  a .  ft  <  I 

It  is  evident  that  equation  (7)  has  faster  reaching  speed  than  that  of  SMC  since  an  exponential  term 
exists  in  equation  (7).  Apart  from  that,  when  |S|  grows  larger,  the  coefficient  of  the  reaching  law  approaches 

K  which  is  larger  than  K  ,  hence,  the  reaching  speed  is  faster.  Meanwhile,  when  the  trajectory  reaches  near 

N(S)  P 

the  sliding  manifold  and  \S\  reaches  0,  then  the  chattering  phenomenon  is  reduced  due  to  the  existence  of  the 
exponential  term  and  the  increase  of  N{s)  . 


3.  DFIG  WIND  TURBINE  APPLICATION 
3.1.  Modeling  of  the  wind  turbine  and  gearbox 

Consider  a  wind  turbine  provided  with  blades  of  length  R  resulting  in  a  generator  through  a 
multiplier  of  speed  of  gain  G 

The  aerodynamic  power  appearing  at  the  level  of  the  rotor  of  the  turbine.  It  is  given  by  [18]: 

P,  =  X-Cp(X).p.n.R1V^  (9) 

where  p  is  the  air  density,  R  is  the  blade  length  and  V  is  the  wind  velocity.  Knowing  the  speed  of  the 
turbine,  the  aerodynamic  torque  is  this  directly  determined  by: 

Taer  =  Jf  (10) 


The  multiplier  of  speed,  which  is  arranged  between  the  wind  turbine  and  the  generator,  aims  at 
adapting  the  speed  of  the  turbine,  rather  slow,  to  the  one  that  requires  the  generator.  It  is  modelled  by 
following  tows  equations 


Q.t 


G 


(ID 
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where  T  is  the  driving  torque  of  the  generator  and  0.mec  is  the  generator  shaft  speed. 

The  power  coefficient  represents  the  aerodynamic  efficiency  on  the  wind  turbine.  It  depends  on  the 
characteristic  of  the  turbine  [19], [20].  The  ratio  of  speed  is  defined  as  the  report  between  the  linear  speed  of 
blades  and  the  wind  speed: 


2  = 


Qt.R 

V 


(12) 


A  typical  relationship  between  C  p  and  A  is  shown  in  Figure  2.  It  is  clear  from  this  figure  that  there 
is  a  value  of  A  for  which  C p  is  maximum  and  that  maximizes  the  power  for  a  given  wind  speed.  The  peak 
power  for  each  wind  speed  occurs  at  the  point  where  C p  is  maximized  [20],[21], 


Figure  2.  Power  coefficient  for  the  wind  turbine  model 


3.2.  Modeling  of  the  Double  Fed  Induction  Generator 

The  model  of  state  of  Park  of  the  generator  under  the  following  matrix  shape  [21]: 


vsd 


Vfd 


v 


rq 


~  ^shd  fisd  ~  ^s&sq 

~  ^shq  fisq  ^s&sd 

—  Rrird  +  —^rfisq 

=  Rrirq  +  —  (j)rq  +0)r(j>rci 


(13) 


where  f?sand  Rr  are  the  stator  and  rotor  phase  resistances,  respectively.  cos  =co+COr ,  and  &>= /?-£lmefis  the 
electrical  speed  where  p  is  the  pair  pole  number.  The  stator  and  rotor  flux  can  be  expressed  as 


\fisd  ~Lsisd  +Mird 
\fisq  =LsiSq  +Mirq 

l^rd  =  Lrird  +  Misd 

\0rd  =  Lrirq  +Misq 


(14) 


(15) 


where  isd ,  i sq ,  ird  and  irq  are  the  direct  and  quadrature  stator  and  rotor  currents,  respectively.  The  active  and 

reactive  powers  exchanged  with  the  electrical  network  also  depend  on  electric  variables  of  the  stator  and  on 
the  rotor.  In  the  stator,  the  active  and  reactive  powers  are  respectively  given  by  [21]: 
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j PS  vsd^sd  vsqfq 
\  Qs  =  vsqhd  ~vsd^sq 
The  electromagnetic  torque  is  expressed  as 

Pem  ~  Pifisdhq  (frsqisd) 


(16) 


(17) 


4.  CONTROL  STARTEGY  OF  THE  DOUBLE  FED  INDUCTION  GENERATOR 
4.1.  Decoupling  of  the  active  and  reactive  powers 

The  principle  of  orientation  of  flux  statorique  consists  in  aligning  the  statorique  flux  according  to 
the  axis  "  d "  by  the  rotating  reference  table.  This  last  constraint  is  favorable  to  have  a  simplified  model  of 
command.  By  setting  the  quadrate  component  of  the  stator  to  the  null  value  as  follows: 

(j)s  —  (ps(i  =>  (f)Sq  —  0  (18) 

Equation  (6)  is  simplified  as  indicated  below: 

\0sd  =Lsisd  +Mird 

[  0  —  LsiSq  +  Mirq 


By  neglecting  the  stator  resistance  Rs 

\vsd  =0 

\vsq  =Vs 


(20) 


Using  equations  (12),  (13),  (14)  and  (17)  the  stator  active  and  reactive  power  can  then  be  expressed 
only  versus  these  rotor  currents  as: 


M 

P  =V  i  =  -V  —  i 

1  s  r  slsq  s  L 


Q,  =v,itJ 


VSM 


lrd 


(21) 


4.2.  SMC  Regulator  Synthesis 

The  Figure  3  represents  the  plan  of  the  simulation.  The  machine  is  directly  connected  to  the  network 
of  the  side  statorique.  Of  an  other,  the  circuit  rotorique  is  fed  through  a  converter  to  MLI.  The  errors  between 
the  reference  powers  and  those  measured  at  the  level  of  the  stator  are  handled  by  the  algorithm  of  considered 
control,  to  design  the  reference  tensions  of  the  rotor.  These  reference  tensions  as  well  as  those  in  the  entrance 
of  the  converter  to  MLI  are  used  by  the  technique  of  modulation  considered  for  the  synthesis  of  the  signals  of 
order  for  the  bidirectional  switches  of  the  converter. 

In  this  section,  we  develop  a  law  command  to  control  the  stator  active  and  reactive  power. 

Choosing  the  following  switching  functions: 

Sps=(Ps-P)j+^\(Ps-Ps)lt 

<  f  \  Y(  \  (22) 

sQs={Qs-Qx)+YjU-Qsk 

where  Ps  ,  Qs  denote  the  reference  stator  active  and  reactive  power,  respectively;  /. p ,  Aq  are  the 
coefficients  of  the  switching  function.  The  derivatives  of  both  switching  functions  are: 
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\sPs  =-ps+iM-p\ 

j  =-Qs  +  'Lq  \Qs  ~  Q  s  ) 


Based  on  equation  (21),  we  have: 

Sps  =  Vs  —  (irq  +  APi  rq )+  APPS 


M  : 


SQS-Vs~T  'rd+AQird-^Q—  +AqQ 


Considering  equation  (13)— ( 14)  and  (15),  we  have 


ird  —  ^  {}’rd  R-Ard  A-  oLrcorirq^ 


:  _  1 
lrq~^Lr 


V rq  Rrbq  tr{  V, 


gM 


(23) 


(24) 


(25) 


where  M 

LsLr 

(24),  we  have 


is  the  scatter  coefficient  of  the  magnetic  flows.  Substituting  equation  (25)  into  equation 


S  = 


Gi 

Pi 


+  [/t 


where 


VSM 
<jLS  Ly 


Gl  -flhd  +  hh-q  +/4 


h=-vs 


MR,- 
oLs  Lr 


h  = 


VSM 

h 


[(Oy 


f  A  ~^pPs 


fn  - 


V*M 

CjL^  l~Jy 


G2  ~  fllh-d  +f33bq  +  f AA 


fll  ~ 


MVS 

4 


(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 

(34) 

(35) 
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/33  “  ~Vs  ' 

/44  =  -V's 


RrM 

gM 2 

oL2Lr 

s 


+\>Qs 


(36) 

(37) 


Now  we  substitute  the  reaching  law  in  equation  (7)  into  equation  (26),  which  yields 
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Figure  3  shows  the  block  diagram  of  the  system  implemented  with  the  sliding  mode  control  (SMC) 
with  exponential  reaching  law  control  (ERL). 


Figure  3.  Configuration  of  the  power  control  of  the  DFIG  with  SMC  with  ERL 


5.  RESULTATS  AND  DISCUSSION 

To  test  the  performance  of  the  proposed  ERL  based  SMC  in  wind  turbine  active/reactive  power 
regulation;  a  comparative  study  is  conducted  under  variable  wind  speed  scenarios.  A  7.5  KW  DFIG  wind 
turbine,  simulation  model  is  constructed  in  Matlab/Simulink.  The  parameters  utilized  are  listed  in  Table  1. 


Table  1.  DFIG  nominal  parameters 

Table  2.  Turbine  parameters 

Table  3.  (DFIG+TURBINE) 

7.5  KW  380/660  V 

50  Hz 

Diameter=  13 

parameters 

P=2  Rs=0.455  Ohm 

Rr=0.62  Ohm 

Gearbox=35 

J=0.3125  kgm2 

M=0.078  H  Ls=0.084  H 

Lr=0.08 1  H 

Number  of  blades=3 

f =0.00673  Kg.m 
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In  this  section,  we  have  chosen  to  compare  the  performances  of  the  DFIG  with  two  different 
controllers.  The  Sliding  Mode  Control  will  be  first  tested  and  will  be  the  reference  compared  to  the  other: 
Sliding  Mode  Control  with  exponential  reaching  law  control  (ERL). 

Figures  4  and  5  show  the  simulations  results  of  the  two  methods  in  electromagnetic  torque  control, 
active  and  reactive  power  tracking  and  current  statorique  and  rotorique.  What  needs  to  be  noted  is  that  the 
parameters  utilized  in  the  programming  of  the  SMC  that  the  Aps  and  Aq  are  both  chosen  to  be  25  and  p  is 

fixed  to  be  0.5,  and  a  is  1  and  the  coefficients  of  the  new  switching  function  in  equation  (22)  are  both 
chosen  to  be  5.  The  solvers  are  still  Runge-Kutta  (ode45)  and  the  sampling  rate  is  25  ps.  As  it  can  be  figured 
out  from  Figures  4  and  5  that  the  performance  of  the  proposed  ERL  based  SMC  is  better  than  that  of  the 
classical  SMC  in  tracking  the  electromagnetic  torque.  Active  and  reactive  power.  Furthermore,  the  chattering 
phenomenon  is  largely  reduced.  It  is  evident  that  the  curve  of  the  ERL  based  SMC  tracks  the  reference  value 
precisely  with  little  chatters.  When  utilizing  the  conventional  SMC  approache,  high  frequency  chatters  with 
high  amplitude  exist.  These  chatters  are  caused  by  the  discontinuous  switching  function  of  the  control  input. 

When  the  ERL  based  SMC  approach  is  utilized,  the  discontinuous  components  of  the  control  input 
have  relatively  higher  gains,  making  the  chatters  more  distinctive.  The  discontinuous  component  of  the 
control  input  automatically  decreases  to  an  acceptable  level  and  the  sliding  stage  becomes  much  more 
smooth  and  precise. 
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Figure.  4  power  control  of  the  DFIG  with  Figure.  5  power  control  of  the  DFIG  with 

HOSTSM  F  -HOSTSM 

6.  CONCLUSION 

In  this  paper,  we  have  presented  a  complete  system  to  produce  electrical  energy  with  a  doubly-fed 
induction  generator  by  the  way  of  a  wind  turbine. 

Once  the  complete  model  of  the  wind  chain  is  validated,  we  spent  optimizing  active  and  reactive 
power  injected  to  the  grid.  In  this  context,  two  robust  control  strategies  have  been  studied  and  compared  to 
show  the  improvement  of  the  quality  of  energy  and  energy  efficiency,  particularly  under  the  influence  of  the 
variation  of  wind  parameters  of  the  wind  turbine. 

The  control  of  the  machine  has  been  presented  in  order  to  regulate  the  active  and  reactive  powers 
exchanged  between  the  machine  and  the  grid.  We  proposes  a  DFIG  wind  turbine  control  technique  with  a 
novel  sliding  mode  control  technique  under  variable  wind  speed  based  on  novel  exponential  reaching  law 
(ERL).  The  proposed  SMC  is  supposed  to  reduce  the  level  of  chattering  phenomenon  and  accelerate  the 
reaching  process,  which  has  been  verified  by  the  numerical  simulation  of  a  nonlinear  controlled  system. 
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